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The QCD phase diagram is studied in the lattice QCD simulation with the imaginary chemical potential 
approach. We employ a clover-improved Wilson fermion action of two-flavors and a renormalization-group 
improved gauge action, and perform the simulation at an intermediate quark mass on a 8 3 x 4 lattice. The QCD 
phase diagram in the imaginary chemical potential fii region is investigated by performing the simulation for 
more than 150 points on the (/3, fii) plane. We find that the Roberge- Weiss phase transition at fii /T — n/3 is 
first order and its endpoint is second order, which are identified by the phase of the Polyakov loop. We determine 
the pseudo-critical line from the susceptibility of the Polyakov loop modulus. We find a clear deviation from a 
linear dependence of the pseudo-critical line on fjij . 
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I. INTRODUCTION 

The QCD phase diagram, which includes states of 
matter formed in terms of the strong interaction, has 
been of prime interest in recent physics covering par- 
ticle physics, hadron/nuclear physics and astrophysics. 
Because QCD is non-perturbative in most regions of 
the QCD phase diagram, one is forced to use the lattice 
QCD in order to obtain a quantitative understanding. 
The lattice QCD is expected to provide reliable infor- 
mation on the phase structure based on QCD. Indeed, 
recently, there have been many active quantitative in- 
vestigations about the finite temperature QCD 

On the other hand, simulations of systems with non- 
zero quark chemical potential /i have been a long chal- 
lenge for the lattice QCD because of the notorious 
sign problem. In the lattice QCD, a fermionic deter- 
minant det A(/i) is used as a probability in a Monte 
Carlo method. The introduction of non-zero [i makes 
det A(/i) complex, and therefore leads to the break- 
down of the stochastic part of the lattice QCD, see 
Ref. H. 

Despite of the severe sign problem, several ap- 
proaches have been proposed to study the QCD with 



nonzero fi, see e.g. [3, 4]. One idea is to avoid the sign 
problem by performing simulations in systems with an 
imaginary chemical potential. A partition function and 
its free-energy are analytic within one phase even if 
chemical potential is extended to complex, which is 
true until the occurance of a phase transition. This val- 
idates the imaginary chemical potential approach for 
the study of the QCD phase diagram. 

Fermion determinants satisfy a well-known relation 



= 75 A(- A1 *) 7 5, 



(1) 
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which holds for a complex chemical potential: /i = 
1*>B, + ifJ>i(fJ>R, fJ>i £ K). This implies that det A(/i) 
is complex for a real chemical potential fi = fin, 
which causes the sign problem. On the other hand, one 
can easily prove det A(/i) is real for a pure imaginary 
chemical potential /i = i/ij. The sign problem does not 
occur in this case, and Monte Carlo methods are avail- 
able. The imaginary chemical potential approach pro- 
vides an insight into the QCD phase diagram through 
the analytic continuation. In addition, data obtained 
in such a simulation can be used for the matching of 
phenomenological models such as Polyakov loop ex- 
tended Nambu-Jona-Lasinio(PNJL) models with the 
lattice QCD flU. 

The imaginary chemical potential approach has been 
studied by using staggered fermions with two flavor 
Igl], three flavor tfToh . four flavor 111 -14 1 in 2-color 
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FIG. 1: Schematic figures for the Nf = 2 QCD phase diagram in the (fi 2 , T) plane (left) and 
/T, T) plane (right). A : Pseudo-critical point at fi = 0. B : Critical endpoint. C : Roberge- 
Weiss endpoint. AB : Pseudo-critical line. AC : Extension of the line AB into the imaginary 
chemical potential plane. CD : Roberge-Weiss phase transition line hi /T = n/3. In the right 
panel, larger /T region of the phase diagram is obtained from the RW periodicity. 



16], by using Wil- 



QCD and finite isospin QCD 
son fermions with two flavor 

Staggered fermions of the standard type might have 
suffered from two problems. First, it needs a fourth- 
root trick for one flavor 1I30I . Second, it does not show a 
scaling behavior expected from three-dimensional Q(4 ) 
spin models for the finite temperature transition IU8II . 
although the possibility of the first oder phase transi- 
tion for two degenerate flavors was studied in Ref. 111911 . 
Wilson fermions are free from the fourth root tricks 
and show the correct scaling behavior. On the other 
hand, Wilson fermions suffer from an explicit break- 
ing of chiral symmetry. However, one can define a 
subtracted chiral condensate, which satisfy a correct 
scaling behavior. Although Wilson fermions require 
more computational time than that required in stag- 
gered fermions, simulations with Wilson fermions are 
now possible even on the physical quark masses at zero 
density. 

In finite temperature simulations with the combina- 
tion of the plaquette gauge action and the standard Wil- 
son quark action at Nt = 4, the transition is smooth 
crossover at small and large quark mass and rapid 
crossover at intermediate quark mass [20), which is 
different from what is expected the transition becomes 
sharp for light and heavy quark masses. This unex- 
pected behavior is removed by improving the gauge 
action 1181 12 ill . The improvement of the gauge action 
is essential in removing lattice artifacts at finite lattice 
spacings. 



Thus, the study of the Wilson fermions with im- 
proved terms is complementary and useful to confirm 
results obtained in other actions l!7Ml3[ Il7ll and to es- 
tablish a better understanding of the QCD phase di- 
agram. In this paper, we study the two-flavor QCD 
phase diagram at an intermediate quark mass by us- 
ing the imaginary chemical potential approach. This is 
the first employment of the two-flavor Wilson fermion 
with a clover term and the renormalization-group(RG) 
improved gauge action to the imaginary chemical po- 
tential approach. 

This paper is organized as follows. In the next sec- 
tion, we briefly review properties and issues of the 
imaginary chemical potential region of the phase di- 
agram. The setup for the simulation is also explained 
here. The numerical results are shown in Sec. |TTT] We 
investigate the deconfinement transition and Roberge- 
Weiss endpoint in Sec. lIII Al and Roberge-Weiss phase 
transition line in Sec. IIIIBI We determine the pseudo- 
critical line in Sec. lIII CI The final section is devoted to 
a summary. 

II. FRAMEWORK 
A. Phase diagram with imaginary chemical potential 

We begin with a brief overview of the QCD phase 
diagram and of issues in question. 

First we show in the left panel of Fig. Q] an expected 
phase diagram in (/i 2 , T) plane containing both the real 
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(a* 2 > 0, fi = fin) and imaginary (/i 2 < 0, fj, = 
ifii) regions. Even if /i 2 < 0, it is expected quark- 
gluon-plasma(QGP) and hadronic phases exist at high 
and low temperatures, respectively. The two phases 
are separated by the deconfinement crossover transi- 
tion line, which is an extension from the /i 2 > re- 
gion. These are consequences of an analyticity of the 
grand partition function. The absolute value of the 
Polyakov loop is often employed to identify confine- 
ment/deconfinement phase, although it is not a real or- 
der parameter because of the crossover nature of the 
transition. 

Two characteristics of the /i 2 < region are 
so-called Roberge-Weiss(RW) pha se transition and 
Roberge-Weiss(RW) periodicity 02211 . The QCD grand 
partition function has a periodicity with a period 

2tt/N c as 



T N c ) 



(2) 



where k is an integer. Furthermore, Roberge and Weiss 
showed from a perturbative analysis the existence of a 
first-order phase transition on the line ////T = ir/N c , 
and from a strong coupling analysis the absence of 
such a transition at low temperatures. These features 
hold for SU^A^) gauge theories. Hereafter we con- 
sider N c = 3. The RW phase transition relates to the 
Z(2>) symmetry and an order parameter identifying this 
phase transition is the phase or imaginary part of the 
Polyakov loop. Because the RW phase transition oc- 
curs at high temperatures but does not at low temper- 
atures, it may have an endpoint at a temperature Trw 
on the line fii/T = tt/3. These features are well man- 
ifested in the (fij/T, T) -phase diagram, see the right 
panel of Fig. [T] 

Although the phase diagrams in Fig. Q] are naively 
expected, several points remain as issues which should 
be discussed further. Figure[T]is drawn according to the 
two points : 

• The RW phase transition exists at high tempera- 
tures and has an endpoint. 

• The extension of the crossover line exists in the 
/i 2 < region. 

The second point needs an assumption that a pseudo- 
critical line can be defined for a crossover. Once the 
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FIG. 2: Counter examples against the naive 
expectation. R is a convergence radius of the 
pseudo-critical line determined in the fj 2 < 
region. 



pseudo-critical line is defined, the second point is en- 
sured by the identity theorem. 

Based on the above two points, one can consider 
several counter examples, against the naive expectation 
Fig.G] see Fig. [2] In order to reveal the phase structure 
without assumptions on imagination, we need to know 

• the location and order of the RW endpoint 

• the location of phase boundary of the deconfine- 
ment transition 

• the way with which the RW and deconfinement 
transition lines are connected 

Indeed, the quark-mass dependence of the RW end- 



s dependence 
. DT 



8|] a possibility that other first-order phase transition 
lines depart from the RW endpoint, which corresponds 
to the left top panel in Fig [2] 

The above questions are important for two reasons. 
First, they are relevant with the definition range of the 
pseudo-critical line determined by the imaginary chem- 
ical potential approach, which also relates to the appli- 
cable range of the line in the /i 2 > region. For in- 
stance, if the pseudo-critical line has an endpoint near 
T pc at jtt = J3TI1 . then a convergence radius R of the 
pseudo-critical line obtained from the imaginary chem- 
ical potential approach is given by a distance between 
the endpoint and T pc . The pseudo-critical line obtained 
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can be applied to a domain with the convergence ra- 
dius. 

Second, it is speculated from the Lee- Yang theo- 
rem lEHE^ that a phase transition in the < region 
is relevant to one in the /i 2 > region, although the 
distribution of the Lee-Yang zeros of the QCD grand 
partition function has not been well understood. As- 
suming the Lee-Yang zeros of the QCD are distributed 
on a line in the complex fugacity plane, it is possible 
that the RW phase transition line and its endpoint in 
the /1 2 < region reflect the first order phase transition 
line and critical endpoint in the /1 2 > region of the 



QCD phase diagram. 



B. Formulation and setup 



We employ the RG-improved gauge action 02511 



S -I 



X,fA<,L> X,fA,l/ 



,(3) 



with ci = —0.331 and Co = 1 — 8ci, and the clover- 
improved Wilson fermion action with the quark matrix 



A(x,y) = 6 X , X , - [(! -7i)Ui(x)5 x>tX+ t + (l + li)Uj(x')S x , !X _ } 

i=i 

- « [e +/1 (l - 74)^)^+4 + e -"(l + 74)^(^)^,»-4 

- cCsif^,!' £ ViivF^v (4) 

I 



Here /1 is the quark chemical potential in lattice unit, 
which is introduced to the temporal part of link vari- 
ables. 

In order to scan the phase diagram, simulations were 
done for more than 150 points on the (/ij, f3) plane in 
the domain < m < 0.28800 and 1.79 < /3 < 2.0. 
Note that the RW phase transition line in the present 
setup is given by /ij = tt/12 ~ 0.2618. All the simu- 
lations were performed on a iV 3 x N t = 8 3 x 4 lattice. 
The value of the hopping parameter k were determined 
for each value of j3 according to a line of the constant 
physics with mps/my = 0.8 obtained in Ref. 
The coefficient of the clover term Csw was determined 
by using a result obtained in the one-loop perturbation 
theory [27] : C S w = (1 - 0.8412/T i)" 3 / 4 . 

The hybrid Monte Carlo algorithm were employed 
to generate gauge configurations. The setup for the 
molecular dynamics was as follows: a step size St = 
0.02, number of the molecular dynamics N T = 50 and 
length N t 6t = 1. The acceptance ratio for this setup 
was more than 90 %. We generated 11, 000 trajecto- 
ries for most parameter sets, and 16, 000 trajectories 
for some parameter sets near the deconfinement tran- 



sition at fj,j = 0. For all the ensemble, the first 5,000 
trajectories were removed as thermalization. The pla- 
quette P, Polyakov loop Le 1 ^ and their susceptibilities 
were measured for each trajectory, where L and <f) are 
the modulus and phase of the Polyakov loop. 

The Polyakov loop operator is as usual defined by 

V c X t=l 

where Ny = iV 3 . The modulus and phase are, after 
the ensemble average, defined by (P i) = Le 1 ^. Their 
susceptibilities are also defined by 

XL = N V ((L-(L))) 2 , (6) 
X^AV((0-W)} 2 . (7) 

Distributions of the Polyakov loop in the complex 
plane are shown in Fig. [3] As the figures clearly show, 
the phase structure can be identified by considering the 
j3 and /ij dependence of the Polyakov loop. 
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FIG. 3: Scatter plots of the Polyakov loop. Left : f3 = 1.80 (low temperature (below T pc )). 
Right : /3 = 1.95 (high temperature (above Tew))- 
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FIG. 4: The /3-dependence of the plaquette P (left) and its susceptibility \p (right) f° r various [ii. 



III. NUMERICAL RESULTS 
A. Deconflnement transition and RW endpoint 

First, we investigate the deconfinement transition 
and RW endpoint, by considering the /3-dependence of 
the observables. 

The plaquette P and its susceptibility xp are shown 
in Fig.@] P is a smooth increasing function of (3 for all 
/ij. The effect of /ij suppresses P for intermediate f3, 
while does not change it for small and large f3. How- 
ever, the effect is up to a few percent. xp has a broad 
peak, and the peak position moves toward larger f3 with 
increasing /i/. 

The Polyakov loop modulus L and its susceptibility 
Xl are shown in Fig. [5] L increases slowly for small (i. 
At a certain value of /?, the slope of L becomes large, 
and xl has a broad peak. These behaviors suggest 



the possibility that the system undergoes the crossover 
transition with increasing f3 or temperature. The peak 
position tends to move towards a larger f3 with increas- 
ing /!/, similar to the behavior of \P- This behav- 
ior confirms that pseudo-critical temperatures become 
higher with the increase of [ij until fij = it/ 12. 

It should be noted that the crossover behavior is ob- 
served for all hi, even on the line /i/ = 7r/12. Hence, 
the pseudo-critical line starts from fii = and reaches 
IH = 7T/12. 

Also note that the peak of xp and xl are not sharp, 
and the signal is unclear. This may come from the 
small spatial and temporal sizes. In fact, the WHOT 
collaboration reported in a finite temperature simula- 
tion with the same action that the hopping parameter 
dependence of the Polyakov loop susceptibility shows 
a pronounced peak in a 16 3 x 6 lattice [26]. 

The Polyakov loop phase <fi and its susceptibility x<t> 
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are shown in Fig. [6] <f> rapidly changes near j3 = 1.92 
only for /ij = tt/12, while it is a smooth function of j3 
for < [ii < 7r/12. It is seen that for /ij = tt/12 there 
is one vacuum at low temperatures and are two vacua at 
high temperatures. The histogram of </> at /i/ = 7r/12in 
Fig. [7] also shows this behavior. The susceptibility x<p 
shows a divergent-like behavior near j3 = 1.92 only 
for fii = tt/12. These behaviors suggest the possibil- 
ity that the system undergoes the second order phase 
transition at the RW endpoint with increasing tempera- 
ture. 

The transition point of x<j> m Fig. [6] agrees with the 
the peak position of xl in Fig. within error bars. 
Hence, the pseudo-critical line is connected with the 
Roberge-Weiss phase transition line at the RW end- 
point. This is the case shown in Fig. [T] We observed 
that near the RW endpoint, the Polyakov loop modulus 



shows the crossover-like behavior and its phase shows 
the second order-like behavior. Note that the order of 
the phase transition were naively obtained from the be- 
haviors of the observables. In order to confirm the or- 
der of the phase transition, the finite volume scaling 
analysis should be investigated in future works. 

B. RW phase transition 

Next, we show the /^-dependence of the observables 
in Figs. [8] and [9] and investigate the nature of the RW 
phase transition line. Results of P are not shown there, 
because the effect of [ij changes the value of P up to a 
few percent, as shown in the previous subsection. 

L decreases for [ij < tt/12, and increases for fij > 
tt/12. A rapidly decreasing behavior, which occurs due 
to the intersection with the pseudo-critical line, is ob- 
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FIG. 7: The histogram of at /xj = tt/12 for 
various /3. 



the RW endpoint depends on the quark mass, and first 
order for light and heavy quark masses and second or- 
der for intermediate quark masses. Hence the second 
order nature of the RW endpoint comes from the inter- 
mediated quark mass. 



C. Pseudo-critical line 



As we have discussed in Sec. IIII Al the Polyakov 
loop modulus L shows the deconfinement crossover 
with increasing temperature, which spans the range 
< [ii < 7r/12. We extract the value of (3 pc by fitting 
five or six data of xl near the peak with a Gaussian 
function: xl °c exp(-6(/3 - (3 pc ) 2 ). 
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TABLE I: The values of /3 pc determined by 
fitting xl with a Gaussian function. 
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FIG. 8: The ^/-dependence of L for various ft. 



served for /3 = 1.87 and 1.90. The line for /3 = 1.95 
crosses the RW phase transition line at n\ = tt/12. 
However, a critical behavior is not observed there. 

The iij -dependence of <f) and x<t> are shown in Fig. [9] 
<j> is a smooth function of [ij at low temperatures (fi = 
1.80-1.90), while jumps to — 2tt/3 from at = 
tt/12 at a high temperature (J3 = 1.95). The system 
undergoes the first order phase transition at fj,i = tt/12 
at high temperatures. Note L is periodic and cj) is anti- 
periodic, which is caused by the periodicity of the 
dependence of the Polyakov loop fl]. 

Together with the result obtained in the previous 
subsection, we find that the RW phase transition is 
the first order one and ends at second order endpoint, 
which are identified by the phase of the Polyakov loop. 
This feature is consistent with the results obtained in 
Ref. [ntl . It was reported Jgl [lot that the order of 



The results are shown in Table U and plotted as a 
function of fij in Fig. [10] The data points in Fig. [10] 
clearly deviate from a linear function near [ij ~ tt/12, 
which implies the contributions of higher-order terms. 
The need of terms of order higher than fi 2 in the func- 
tional form of the pseudo-critical line was first pointed 
out in ifnUlrjll . Such a behavior has not been obtained 
from a study with the combination of plaquette gauge 
action and the standard Wilson fermion without the 
clover term J17I1 . Hence the clover- improved Wilson 
fermion and the RG-improved gauge action leads to 
the deviation from a linear dependence of the pseudo- 
critical line on [ij. This finding is an advantage of the 
improved actions. 

In general, the pseudo-critical line can be expanded 
in powers of /ij : 



2\n 
I) ) 



(8) 
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FIG. 9: The jii -dependence of <f) and \<P f° r various f3. 
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FIG. 10: /3 C and /ij. 



which is defined for the range < fii < w/ 12 bounded 
by the RW endpoint. Within the definition range and 
with the present numerical results, only a few terms 
can be determined. Here, we test quadratic and quartic 
functions and a Pade approximation of a simple type 



1 



Ppciin) = coy 



(9) 



The results are shown in Table UT1 and Fig. QT| The 
Pade approximation and quartic function are better 
than the quadratic function because of non-/i| contri- 
butions, as we have mentioned above. The Pade ap- 
proximation is slightly better than the quartic function 
in x 2 /d.o.f, but the difference is quite small in the 
fi 2 < region. They are extended to the /i 2 > plane 
through the analytic continuation /ij — > —fi 2 . The re- 
sult in the fi 2 > region depends on the fit functions 
in particular at large /i 2 . 

The pseudo-critical line j3 pc (fii) can be transformed 
into the one in physical unit through the relation 

1 



T = 



(10) 



aQ9)JV ( ' 

However, we use the data for (3 dependence of 
T/T pc (0) obtained in Ref. M26I1 instead of determining 
lattice spacings. Note that the value of f3 pc (0) slightly 
disagrees with the one obtained in Ref. 12611 . which 
causes the deviation of T pc /T pc (0) at \i = from 1 
about 1% 

TABLE III: The coefficients of the fit functions for the 
pseudo-critical line. 



TABLE II: The coefficients of the fit functions for the pseudo- 
critical line. Co is fixed with the central value of /3 pc (0). 



type c 


Cl 


C2 


Vx 2 /d.o.f 


quadratic 1.866 


0.67(4) 




1.56 


quartic 1.866 


0.40(9) 


6.64(2.17) 


0.96 


Pade 1.866 


-6.59(90) 


-6.84(88) 


0.76 



type dp di da <fe \/x 2 / d -°-.f 

quadratic 1.01 0.077(5) - - 1.68 

quartic 1.01 0.039(12) 0.060(19) - 1.00 

Pade (I) 1.01 -0.44(5) -0.49(5) 0.70 

Pade (II) 1.01 -0.89(10) -1.02(12) 0.119(37) 0.55 

It should be noted that T is not a linear function of (3, 
therefore there is no need that the above functions can 
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FIG. 1 1 : The pseudo-critical line f3 pc in the imaginary(left panel) and real(right panel) region. 




FIG. 12: The pseudo-critical line in dimensionless physical unit in the imaginary (left panel) 
and real(right panel) region. 



be used here. We consider quadratic, quartic functions 
and two types of the Pade approximation. 

2n 

(ID 



TO ~ a " 

P c n 



Ml 



T, 



pc 



^=4 , , rT7 / \ 9 (Pade (I)), (12) 



^pc 



pc 
2 



= rf 



l + d 2 (/i//T pc ) 2 

1 + rfi(/t//T 3 



l + d 2 (A//T pc ) 2 + d 3 (/i//T P c) 4 
(Pade (II)), (13) 

where /ij = a/ij, and /}/ is the imaginary chemical po- 
tential in physical unit. T° c and T pc are pseudo-critical 
temperatures at zero and finite chemical potentials. As 
we have mentioned above, drj(= T pc /T pc (0) at fi = 0) 
deviates from one with 1 % because of the disagreement 
of f3 pc (0) from Ref. Q. Here we added Pade (II) de- 
fined in Eq. ( fT3l , which was investigated to fit the criti- 
cal line for four flavors in Ref. [14]. Note that Eq. (flj) 



can be extended to T pc = for the /i 2 > region II14I1 . 
although the definition range of the line is restricted by 
the RW endpoint. The results are shown in Table Hill 
and Fig. fT2l 

The central values of the quadratic, quartic and Pade 
(I) approximation show similar behavior to /3j, c (/•*/)■ 
The Pade (II) are consistent with the quartic and Pade 
(I) until Hi /T pc < 0.8, and shows sharp rising for 
0.8 < Hi/T pc , Considering the errors, the quadratic 
function still undershoots the obtained data, and the 
quartic function suffers from large errors. The two 
Pade approximations reproduce the data well with 
small errors. The difference between the two Pade ap- 
proximations are observed for 0.8 < Hi/T. Further in- 
vestigations of this region would be important for better 
determination of the pseudo-critical line. 

Similar to f3 pc , the deviation becomes larger with in- 
creasing fi/T in the fj 2 > region and amounts to 
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more than 10% at fi/T ~ 1. The quartic function in- 
creases, and the other three functions decreases. The 
quartic function is completely different from the other 
three functions, although the quartic and Pade (I) are 
almost same in the n 2 < region. 

The quadratic function decreases the fastest, the 
Pade (II) does the next, and the Pade (I) overshoots the 
other two. The curvature at fi/T pc = of a power se- 
ries of (jl/TrT pc ) 2 is often used to make comparisons 
with various studies |28]. We obtain 



t 2 = 7r 2 di = 0.38(12). 



(14) 



Here we employ the result for the quartic function. The 
result is smaller than the value obtained from staggered 
fermions Q7l|29|] and standard Wilson fermions 1U7II and 
implies the pseudo-critical line decreases slower. 

IV. SUMMARY AND OUTLOOK 



RW endpoint and the crossover for the deconfinement. 
The corresponding order parameter for the phase tran- 
sitions are the phase of the Polyakov loop for the RW 
phase transition and RW endpoint and the modulus of 
the Polyakov loop for the pseudo-critical line. The 
Polyakov loop modulus did not show a critical behav- 
ior on the RW phase transition line, while the Polyakov 
loop phase did not on the pseudo-critical line. We de- 
termined the pseudo-critical line from the susceptibility 
of the Polyakov loop modulus. We found a clear devi- 
ation from a linear dependence of the pseudo-critical 
line on fi 2 . 

The present calculation was performed with the in- 
termediate quark mass and small lattice. The finite vol- 
ume scaling analysis and quark mass-dependence anal- 
ysis are necessary to confirm the present results. In par- 
ticular, the order of the RW endpoint depends on the 
mass of the quark. The improvement on these points 
should be done in a future study. 



We have investigated the two-flavor QCD phase di- 
agram in the lattice QCD. The imaginary chemical po- 
tential approach was employed in order to avoid the 
sign problem. The clover-improved Wilson action and 
renormalization-group improved gauge action was first 
applied to the imaginary chemical potential approach. 
The simulation was performed on the 8 3 x 4 lattice 
and at the intermediate quark mass. The simulation 
was performed for more than 150 points in the param- 
eter (/3, fij) plane. Considering the Polyakov loop, the 
imaginary chemical potential region of the phase dia- 
gram was examined. 

Obtained behaviors of the phase transitions are first 
order for the RW phase transition, second order for the 
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